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O ! I. INTRODUCTION 

(N 

In classical Newton-Maxwell physics the electromagnetic (EM) fields have no influence upon gravity, which is 
generated by sources of mass. In general relativity EM fields alter the metric of spacetime and induce a gravitational 
■ force through their energy-momentum tensor 
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is the electromagnetic tensor and is the four-potential. enters the r.h.s. of the Einstein equations 

i?t = kT^. (3) 



O We have taken into account that T'' = 0. The Einstein constant is 

cr. 

5-H . K= —r- = 2.07 X IQ'^^s^ Icm.q, (4) 

where G = 6.674 x lQ~^cw? / g.s^ is the Newton constant and c — 2.998 x 10^°cto/s is the speed of light. We shall 

' use the Gauss system (CGS) of nonrelativistic units and occasionally the international system of practical units. 
, In addition, the Maxwell equations are coupled to gravity through the covariant derivatives of F^y 

^ ' 14^ 



c 

Here g is the metric's determinant, usual derivatives are denoted by subscripts, is the four-current, = dx^ /ds is 
the four-velocity of the charged particles with charge density a. We shall study mainly electrovacuum solutions with 
cr ^ only on some surface, specifying the boundary conditions. The Einstein-Maxwell equations (3,5) show how the 
EM-field leaves its imprint on the metric, which has to satisfy the Rainich conditions Q, 
The gravitational force acting on a test particle is represented by the four- acceleration 



'^p are the Christoffel symbols. At rest vP — (500)^^^^ and 

5p = y(lngoo)^ (7) 
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for an arbitrary metric ^q^. 

In this paper we investigate the problem whether EM-fields can induce strong enough acceleration, rising above the 
gravimeter's threshold of lQ~^cm/s^ or even comparable to the mean Earth acceleration — 980.665cm/s^. Eqs. 
(3,4) show that the metric will be very near to the flat one without any singularities and faraway from the metric 
of a black hole. The question is whether the 20 orders of magnitude supplied by in Eq. (7) are enough to 
lift the EM-gravitational force to that of the Newtonian gravity of very massive bodies. In fact, we should consider 
the contravariant physical (tetrad) four-vector g^^^ = ?7'"'5(,y), where ry^^ — diag (1, — 1, — 1, — 1), but in cartesian 
coordinates and for such an almost flat metric it is indistinguishable from or g(^^-j except for a sign change. 

It seems natural to use perturbation theory in the harmonic gauge where 

gp,!^ = Vi^v + h^^, A/i^^ = -2kTu^ (8) 

and the Poisson equation shows that h^^, is extremely small, while ~ c^k = 1.85 x IQ^^^CGS. 

More concretely, one can ask what is the gravitational acceleration inside a solenoid, the most common device 
for the creation of magnetic fields. There is an exact global solution when the solenoid is infinitely long (^]. It is 
comprised from the Melvin geon as an interior and the vacuum Levi-Civita cylindrical solution Q as an exterior. 
The acceleration has only a radial component and when the metric's signature and the scale of the magnetic field is 
changed correspondingly, reads inside the solenoid 

Hq is the physical magnetic field at the central axis. The acceleration points inwards and reaches its maximum at the 
solenoid's wall r = ri. When ,e.g., ri = Im we have 

5Wmax = 3.69x 10-27i72^ (10) 

Although the dependence on the magnetic field is quadratic, one needs Hq — 1.65 x lO^^G to rise above the gravimeter's 
threshold and Hq = 5.15 x lO^^'G to equal ge- Such strong fields may be present on pulsars, but on Earth the static 
record is 33. 2r (from HFML of the Nijmegen University), while in pulse regime GOT has been obtained. 

Formula (9) has been confirmed by an approximate solution to the Ernst equation for a finite solenoid of spherical 
form 0. Exact solutions were found also for the solenoid's analogues in three and two dimensions and 
expressions similar to Eq. (9) have been derived. It appears that the problem has been solved and the effect is 
negligible. However, there is a loophole, which is discussed in the sequel. 

In Sec. II it is shown that the gravitational acceleration in Weyl-Majumdar-Papapetrou (WMF) fields prLll^lTsLll^ 
includes a term proportional to ^/k and linear in the fields. In Sec. Ill axially-symmetric static Einstein-Maxwell fields 
are revisited. The three classes of Weyl solutions are described, as well as their relation to the general solution and 
situations where they become the most general solution. It is argued that charge and current distributions determine 
their harmonic master-potential and induce them as pure electromagnetic effect upon the spacetime metric. A general 
formula is given for the acceleration and the question of hidden mass sources is investigated. In Sec. IV a short review 
is given of WMP fields in the general static case, in the presence of charged matter, pressureless (dust) or perfect fluid. 
The connection between axially-symmetric static electro (magneto) vacs and stationary vacuum solutions is elucidated. 

In Sec. V we give four electrostatic examples of charged surfaces, which induce Weyl solutions. They involve plane, 
spherical, spheroidal or cylindrical symmetry. The charged plane, the Reissner-Nordstrom (RN) solution and the 
charged cylinder are discussed from a Weyl point of view. The peculiarities in the junction conditions for Weyl fields 
are pointed out. Different physical issues, such as the force arising inside a charged capacitor, repulsive gravity around 
a charged sphere and a chamber with artificial gravity are explored. In Sec. VI we give four magnetostatic examples, 
some of them being analogues to the examples in the previous section. The first two deal with current loops, open 
solenoids and closed spherical or spheroidal solenoids, which induce Weyl fields with enhanced acceleration. The third 
example is the infinitely long solenoid mentioned in the introduction. It has deep relations to Weyl fields. The genesis 
of Eq. (9) is described in detail. The fourth example is the gravitational field of a line current, which is similar to 
Weyl fields with cylindrical symmetry. In Sec. VII the main results about WMF fields are summarised into 19 points. 
The last section contains discussion of their role in physics and some historical remarks. 



II. ROOT GRAVITY 



Let us assume that the metric and the EM-fields do not depend on time. In this stationary case let us further 
simplify the problem by setting = [<}), 0, 0, 0). There is just an electric field 

= Fo^ = (11) 
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Obviously, from Eq. (1) contains only quadratic terms in (j)^. This allows to hide k from Eq. (3) by normalizing 
the electric potential to a dimcnsionless quantity 

/ 1^ _ 

(12) 

The factor Sir is chosen for future convenience. We shall see that this is a much more elegant way to get rid of 
the constants in the Einstein-Maxwell equations than the choice of relativistic units c = 1 and G = 1 (SttG = 1 
sometimes) . 

Imagine now that in some exact solution is proportional to the electric field, contrary to the quadratic dependence 
inEq. (9) 

„2, 




= ac (l)^ = ac J —0^, (13) 



where a is some slowly changing function of order O (1). For example, when a ~ a (.goo) then Eqs. (7,13) lead to the 
functional dependence / = 500 = F (j^)- Let us further assume that the spacetinie is static. Then the above functional 
dependence has the unique form [lll[o.[l3l| 

f=l + Bcf> + ^\ (14) 

In the axially-symmetric case Eq. (14) was found by Weyl already in 1917 0| and such solutions are known as Weyl 
fields. The potential in Eq. (12) is very small everywhere and naturally goes to zero at infinity. Then asymptotic 
flatness fixes the first term which is otherwise an arbitrary constant. It is also fixed by the requirement to go back to 
Minkowski spacetime when (j) — since we are studying the EM-effect on gravity with no masses present. Hence, we 
should not make any gauge transformations (f) (f) + C, which can eliminate the arbitrary constant B but spoil the 
behaviour of the metric. Later we will give arguments that the typical value of _B is 2 and not zero, so that the linear 
term in Eq. (14) is really present. 

Thus in Weyl-Majumdar-Papapetrou fields we have 



3. = ^'r\^^l^^.^rJ<^.)- (15) 

The first term is of the type given in Eq. (13), the second resembles the expression in Eq. (9). Let us note that 




2.58 X 10-^ c^—^^ = 7.37 x IQ-^'. (16) 



Due to the square root, the first coefficient is 10^^ times bigger than the second. We shall call gravitational fields, 
which have acceleration terms ^ -^/k, root gravity. The WMP fields are an example, but there are others too. Thus 
general relativity has a Newtonian limit in the case of mass sources, where ^ G and a Maxwellian limit in the case 
of EM-sources, where ~ ^/G. 

Provided that B k, 2 our search for a strong gravitational acceleration induced by EM-fields doesn't seem so doomed 
as in the introduction. Electrostatic generators can create potential differences of six million volts or 

^^ax - 2 X IQ^GGS. (17) 

If it were applied to a capacitor with distance of \cm between the plates, the electric field will be of the same order. It 
compensates the root coefficient in Eq. (15) and we get acceleration of about lew/s^, \yhich is p erfectly measurable. 
Static magnetic fields create the same gravitational effects as static electric fields Il7j . so a field of 33. 2T 

may induce in principle acceleration of about 10cm/ . One needs just two orders more to counter Earth's gravity. 
These effects are much stronger than any other known general relativistic effects, including gravitational waves and 
gravitomagnetism, which are currently under intensive study. They may be produced in a finite region in a laboratory 
if we learn how to create WMP fields. 

Up to now the only sources of gravitation have been masses (preferably in enormous quantities), therefore, it is 
worth to compare them to root gravity in equal conditions. Let us take a ball of radius R, mass density ji and total 
mass M . The acceleration at its surface is 

9R^^ = = 2.8 X 10- Vi?- (18) 
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Let the substance of the ball be Osmium, which possesses one of the biggest densities, ^ = 22 ASg/cw?. For a 
laboratory set-up R « lO^cm. This yields 6.44 x 10~*cm/s^, which is much less than the root gravity results. 
Comparing Eqs. (14) with B — 2 and (18) we get the effective EM-field necessary to induce the same gravitational 
acceleration 

£;e// = 1.08 X lO-^-R- (19) 

In other words, it is thousand times more effective than a mass source. Root gravity is about 10^^ times stronger than 
usual perturbative gravity arising from Eq. (8), than the exact solution in Eq. (9) and than the second term in Eq. 
(15). Therefore, in the next sections we revisit the WMP fields, putting the emphasis on their physical apphcability. 

III. WEYL FIELDS REVISITED 

Let us start with the axially-symmetric static metric 

ds^ = / {dx°Y - [e^'^ (dr^ + + ^2 ^^2] ^ (20) 

where — ct, — ip, x'^ — r, x^ — z are cylindrical coordinates, / — e^" and u is the first, while k is the second 
gravitational potential. Both of them depend only on r and z. Let — ((^,x,0, O) where x is the true magnetic 
potential. Following Tauber [T^ we introduce the auxiliary potential A 

/ / 

K = -Xz, K = Xr, (21) 

r r 



so that 



p^r^l pzv^J (22) 

r r 

describe the axial and the radial components of the magnetic field. For the electric field one has 

Er = For = — (/ir, — Fq^ = —<j)z- (23) 

The field equations read 

Au = {<I>1 + ^l+Xl + Xl) , (24) 

A(l} = 2{Ur(j)r+U,(l3,), (25) 

AX = 2{UrXr+UzXz) , (26) 

cj)rXz = cj)zXr, (27) 

ul-ul- (02 -cPl + Xl- xl) , (28) 
k 

— = 2UrU, - 2e"^'' {(l}r(t)z + XrXz) , (29) 

r 

krr + kzz = Aw - [ul + ul) , (30) 

where A — drr + dzz + dr/r. We have used the definition given in Eq. (12) and a similar one for A. Using Eq. (27) one 
can prove that A = A ((/)) and the dependence is linear This result holds also for a general static metric [l7l| . 

It is enough to engage just an electric field, there being a trivial magnetovac analogue to every electrovac solution. 
Eqs. (24-27) reduce to ^ 

Au = e'^"" {4>l + 4>l) , Acj) = 2iUr4r+Uz(l)z), (31) 
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which determine 4> and /. Eqs. (28,29) become 

^=ul-ul- (02 _ ^ ^ ^ 2uru, - 2e-2",^^0^ (32) 

and determine k by integration, while Eq.(30) holds identically and is redundant. When = 0, Eq. (31) becomes the 
Laplace equation for u and Eq. (32) gives k (u). These are the axially-symmetric vacuum equations, also discovered 
by Weyl. 

Now let us make the assumption that the gravitational and the electric potential have the same equipotential 
surfaces, / — f {4>)- Eq. (31) yields 

iU^ -2){,pl + ^l)^0 (33) 
and that's how the quadratic relation (14) appears. Replacing Eq. (14) in Eq. (31) one comes to an equation for </) 

We put for definiteness B > 0. Eqs. (14,34) with B < are obtained by changing the sign of 0. The general solution 
of Eq. (34) is not known. However, let us make one more assumption, that depends on r, z through some function 
ip{r,z). Eq. (34) becomes 



{B + 2(b) 0^ A^b 



where j = r, z. Eq. (32) becomes 



(35) 



^=1 l + B^+r - 



<Pi = f^u 0» = /(0)V'», (37) 



^r{'ij)l-^D, = ^riprip^, (38) 



<j>^ 1 + 50 + 02 -02 _,_ ^2 • 

If satisfies the Laplace equation A-0 = 0,0 (-0, B) is determined implicitly T?! from 



A very important equality follows 



where D = B^ — 4. Obviously, k k always, making it much smaller than u. 

Thus in Weyl electrovac solutions the harmonic master potential determines the electric and the gravitational 
fields like u does this in the vacuum case. One may go further and find a relation between and u, transforming Weyl 
electrovacs into Weyl vacuum solutions, although usual transformations work the other way round 0- In particular 
solutions is usually proportional to the charge, = g0. Eq. (36) shows that — qip where is harmonic and finite 
when the electric field is turned off by g ^ 0. In this limit, when B does not depend on q, we have / — > 1 from Eq. 
(14) and fc — > from Eq. (38). Trivial flat spacetime is the result. However, if _B = B/q then / ^ 1 + Bcj), fi B(f)i 
and from Eq. (37) it follows that u = B'iJj/2 and is harmonic. Eq. (38) then gives the vacuum expression for k. 
Hence, we obtain a Weyl vacuum solution with the same (up to a constant) harmonic function ■0. A mass term has 
appeared out of the vanishing charge. 

Let us go back to the electrovac problem. One can add a constant 0o to in order to satisfy the conditions 
■0— >0, / ^ 1,0^0 at infinity or when the electric field is turned off. The integral in Eq. (36) can be analytically 
evaluated and the dependence {ip, B) made explicit. There are three cases, according to the sign of D. The simplest 
one is D = {B = 2). Then / becomes a perfect square and -00 = —1, 

^=-1-^— = -^, f = {l-i^)-\ (39) 

■0 + 00 1-0 

When D <Q Eq. (38) gives —D < 4 and trigonometric functions appear 

2 B 
ipo — —^= arctan , (40) 
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D 



4cos2 {ip + ipo) 




(41) 



When -00 = these formulas coincide with the Bonnor's ones []J|- 

Finally, when D > there exist two expressions for the integral, one as a logarithm, the other in hyperbolic 
functions. They lead to 



coth— (^ + Vo)- —^7 \ ^ , (42) 

2 2 2 B + ^/D-(B- Vd] 




f = = f ^'^^^ 7^ ^i"^ I ' (43) 

4sinh2^(V' + ?Ao) V 

gvTJv-o ^ (44) 

According to Bonnor (who does not introduce ipo) the expressions for / and 4> in the case D > are obtained from 
Eq. (40) by continuation of \/—D to imaginary values and consequently tan — > tanh, cos — > cosh. This, however, 
holds when (2(/) + B) < D , i.e., 4/ < 0, which is unphysical. In the physical case we must also do the replacement 
tanh coth, cosh sinh. The above discussion shows that the point B — 2 has a privileged position, unlike 5 = 0. 

The Weyl solutions were derived with two assumptions imposed on the system: f — f (</)); (j^ — 4' (^)j ^"0 = 0- They 
are particular solutions of Eqs. (31,32). However, when the symmetry is stronger than axial and the fields depend 
on just one coordinate x (not necessarily cylindrical, but a function of r, z), they comprise the general solution. For 
then / (x) and (/) {x) obviously are functionally related and one can always find X [x) so that AX {x) = 0. Taking 

= X {x) and expressing / and </> as functions of X, leads inevitably to the Weyl solutions. In the case of plane 
symmetry x = z, X = x. Cylindrical symmetry gives x = r, X = \n.x. Spherical symmetry has x = vr^^+z^, 
X = l/x. 

We have described the advantages of WMP solutions in inducing a powerful gravitational force. The natural 
questions appear: is it possible to create such solutions in a laboratory? What kind of charged sources should we 
take? Point |20j | and line [U sources lead to singularities and other problems. Therefore, we take a charged closed 
rotationally-symmetric surface with invariant density of the surface charge ■ The electrostatic theorem of Gauss has 
a generalization in general relativity (13i.i22. |. Integrating the r.h.s. of Eq. (5) one obtains the total charge contained 
in some volume 



g£S = / a^cPS, as = aV-gW, (45) 



where is the three-dimensional invariant density and g'^'^^ is the determinant of the space part of the metric. When 
the charge is attached to a surface, one should use tis instead. Integration of the l.h.s. of Eq. (5) leads to a relation 
between e and the electric flux through a closed surface S, encompassing the charged volume 



Aire ~ 



pOl 9{X2,X3) ^ d {X3 , xi ] 



d (u, v) 



d [u, v) 



^ d{xi,X2) 

d {u, v) 



-gdudv. 



For Weyl solutions 



F^^^g = 



(46) 



(47) 



which is the result for flat spacetime. Hence, Eq. (46) becomes the Gauss theorem in classical electrostatics, but with 
(p replaced by i/j, which satisfies the Laplace equation. This fact was already stressed by Bonnor [5lll9|| but in view of 
its importance we have discussed it again. Following a well-known procedure, we obtain a boundary condition on S 
for the jump of the normal component tjjn- 



(48) 
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When ip is given on S, there are two well-defined Dirichlet boundary problems and it may be continued as harmonic 
function inside and outside S. li a < ips < (3, these inequalities hold for throughout space and it will be regular. 
Then /, k and are found from m a manner already explained. The jump of V'n at 5* determines the source Us- 
The inverse is also true. When ct^ is given, there is a unique global t/j, satisfying Eq. (48). Consequently, for any 
distribution of charges on S one can find the electric and gravitational fields they induce. The same can be done when 
S is infinite and / or not closed, but singularities may creep into the solutions. 
Finally, by replacing Eq. (37) into Eq. (15) and expressing f — f we obtain 



For realistic EM-fields / is very close to one and this Maxwellian effect is the only one to be observed. Typical 
Einsteinian effects like light bending, gravitational redshift, time delay, changes in lengths are not enhanced by and 
are negligible. 

Some questions immediately arise. Why do we get a solution for an arbitrary as when the Weyl fields are not 
the most general solutions of Eq. (31)? Why 0, /, k depend not alone on ip but also on the constant B7 How is its 
value determined, can we increase it, to enhance the effect of root gravity? In order to answer them we must return 
to the starting point, Eq.(3). Even when the EM-field is absent, Eq. (3) still has a number of non-trivial vacuum 
solutions, including e.g. gravitational waves. The situation is similar to classical electrodynamics without sources. 
Non-trivial solutions exist, but have to be time-dependent. These are the well-known electromagnetic waves. General 
relativity is a highly non-linear theory and vacuum solutions exist also in the static case. Their sources are well- 
hidden masses and even today there is a gap between the mathematical derivation of solutions [sj and their physical 
interpretation . When EM-fields are turned on, these parasitic masses do not disappear and obscure the pure effect 
of electromagnetism on gravity. Let us try to get rid of them, step by step. First of all, the metric should inherit the 
symmetry of EM-fields. Let us confine again ourselves to axial symmetry. There are a lot of generation techniques, 
which produce non-Weyl solutions of Eq. (31). One of them, "coordinate modelling", adapts the coordinates after 
the equipotential surfaces of the electric field [23 . l25l l2^ and includes in a natural way the set of Weyl solutions. 
Most of the methods (see Ref. Q, Ch.34), however, start from the reformulation of Eq. (31) in terms of the Ernst 
potential E = f — (j)^ [23 , which is real in the absence of rotation. 



The general solution of the Ernst equations can be found when the behaviour of E (z) and / (z) on the axis is given. 
It determines the multipole structure and is useful in astrophysics for modelling the gravitational field of stars. The 
presence of masses is welcomed, since they give the most substantial gravitational effect, followed by rotation (it can 
be incorporated into the formalism), magnetic fields and electric charge at the last place. The Ernst equation is much 
more difficult than the Laplace one and Dirichlet boundary value problems for it were discussed only recently ,28j . 
On the other side, a harmonic function may be easily restored from its values on the axis ip (z) P^|3Q| 




This real expression was used in general relativity for axially-symmetric static vacuum solutions but it is not 
difficult to adapt it to Weyl electrovacs too. In the general solution / (z) is not correlated with (j) (z) and can be 
arbitrary, even when 0(z) vanishes, signalling the presence of masses. In the Weyl solution the metric goes to flat 
Minkowski spacetime when ij: (z) vanishes. Thus the general solution expands over the Weyl one by the addition of 
masses. Probably the same is true for the solutions of Eq. (34), which also contain root gravity terms, since Eq. (14) 
is satisfied. The reason is that Weyl solutions form a complete system, covering the effect of any charge distribution 
and more general solutions can include the only other source of gravitation. More precisely, Weyl fields form an 
overcomplete system due to B which is not fixed by (Tj. In the case of plane, spherical, spheroidal or cylindrical 
symmetry they coincide with the most general solution of the Ernst equation (50). It is quite improbable that the 
unwanted masses should disappear exactly in these cases, so the only way to show their presence is through the value 
of B. Curiously, in a recent paper ^] it is asserted that the sphere, rod and plane are the only Weyl fields with 
geodesic lines of force and are algebraically special of Petrov type D. A logical step is to accept that plays the role 
of (j) in any situation in electrostatics, not only for charged surfaces. Hence, when the gravitation created by charges is 
taken into account, it seems that Weyl fields generalize the solution to classical electrostatic problems. The physical 
electric field is found with the help of Eq. (37) 




(49) 



/Ai; = VfVE, fA(b = V/V0. 



(50) 





(51) 



— — {gooQii) 




'i — 



(52) 
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Now, since /, k are extremely close to one and zero respectively, one can do a perturbation theory around the exact 
Weyl solutions and set Ef^^^ ^ — In the present case WMP fields are similar to instantons, monopoles, solitons and 
other non-perturbative exact solutions in quantum field theory. With high precision all electrostatic formulae hold 
also in the Einstein-Maxwell " already unified theory" . The only new effect is the appearance of an electromagnetically 
induced gravitational acceleration, which reads from Eq. (49), again with high precision, 

= -jVGE,. (53) 

We shall give some arguments in the following that B — 2 when unbiased by parasitic masses. Therefore, as already 
explained, measurable gi are present from the already available electric and magnetic fields. One can reach g,, when 
E, = 1.14 X lO^V/cm = 3.8 x lO^CGS or H, = 380T = 3.8 x lO^G. The fines of acceleration follow the electric field 
lines. Test particles will stay in equilibrium if they are charged and the relation between their mass m and charge e is 



Gm. (54) 

Root gravity has some peculiar features. Changing the direction of Ei one changes the direction of gi and when it 
points upwards with respect to the Earth's surface one has " anti-gravity" . This is true because in our perturbation 
theory accelerations from electric fields and masses like the Earth or laboratory masses are added as usual vectors. 
The exact Weyl solution is necessary to clarify the gravitational induction in a laboratory set-up of finite size in 
space where Ei is present. Although we have used a long-range interaction to induce another long-range interaction, 
in reality static EM-fields are always confined. Eq. (53) shows that putting a Faraday cage on Ei does the same 
on gi and confines root gravity too. It is understandable that when non-mass sources of gravitation are applied, 
the appearance of monopole terms (usually considered as mass terms) is not obligatory. Their existence is usually 
based on the Whittaker's theorem [221, which demonstrates the influence of some combination of the T^j^ components 
(called gravitational mass) upon the gravitational acceleration. The relation is given in terms of surface and volume 
integrals, appearing when the time-time component of the Einstein equations (1) is integrated. In this way it is not 
something separate and additional to them, but a consequence that can't contradict the conclusions following from 
them. Concretely, for Weyl fields this theorem is just the Gauss theorem for grad u. In the case of electromagnetic 
sources the "gravitational mass" is in fact some kind of "energy", inducing gravitational acceleration not necessarily 
with a monopole term. We avoid arguments based on the energy of the gravitational field because its density is not a 
tensor and there are at least five energy-momentum complexes '3Jl|, each with its own merits. Of course, some small 
mass term will always exist, due to the mass of the surface S. However, it will be of usual perturbative nature, many 
orders of magnitude smaller than root gravity. 

Let us turn now to the case of magnetostatics. As was mentioned before, the analogue of (j) is A. It should be 
replaced in Eqs. (14,15,31,37,39,40,42). The analogue of Eq. (46) vanishes because there are no magnetic charges. 
One should take a closed surface with surface current. In an axially-symmetric problem it has just one component, J^p. 
A Weyl magnetostatic solution was given for the first time by Papapetrou in 1947 The analogy with electrostatics 
was investigated by Bonnor ^| who showed that tp is equivalent to the scalar magnetic potential. Then Eq. (48) 
should give the jump of the tangential to the surface component Ht which follows classically from "0 

Ht\t = —J^ (55) 
c 

and is perpendicular to J^p. Eqs. (52,53) still hold with E^i-f H^i-f, so that the gravitational effects of static 
magnetic fields mirror that in electrostatics. One can also imagine a vector potential, corresponding to ^, which is 
more convenient in classical magnetostatics. However, in order to find the metric, the scalar potential tp for each 
classical problem should be calculated too. The fields mainly of linear sources, like a current loop and disks 
|34i l35j| have been examined without noticing the presence of root gravity. 

IV. CONNECTIONS, GENERALIZATIONS AND ANALOGIES 

It is well known that in vacuum spacctimes with isometrics any Killin g ve ctor may be considered as an electro- 
magnetic potential and satisfies the Maxwell equations in this background |36j | . This Papapetrou field is used to find 
EM solutions from vacuum ones ^7|- the axially-symmetric case this leads to dependence between and / and 
one comes to solutions of Eq. (34), which contain ro ot g ravity terms because of Eq. (14). Even a Weyl solution was 
derived recently, starting from the vacuum 7- metric |38| . although this was not stated. 

In the general stationary case the interval reads 

ds^ = f {dx° + Wada;") f-^-/abdx''dx\ (56) 
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where tOa is the gravitomagnetic potential (a = 1,2,3) and ^ab is the three-dimensional metric. In the general static 
case LOa — Q and in its electrostatic subcase the Einstein-Maxwell equations read 

Au = e"2"V0V0, A0 = 2VuV?i, (57) 

where F^a — —4'a- In magnetostatics cj) \ the latter being defined by 

F"'' = (-.9)"^/^e"^'=Ac. (59) 

This equation generalizes Eq. (22). The gradients, the Laplacian and the three-dimensional Ricci tensor are with 
respect to the metric 7ab. Eqs. (57,58) generalize Eqs. (31,32). However, 70;, also enters Eq. (57), making all equations 
interconnected and the system very difficult to deal with. In the special case when Eq. (14) holds with B = 2, ^ab 
becomes flat and Eq. (57) can be solved, since it decouples from Eq. (58). The result is Eq. (39) with a harmonic 
ijj {(p,r, z). Usually one takes 1 — ip = U and = U~^, obtaining the already mentioned Majumdar-Papapetrou 
solutions from 1947 which are conformastatic. Some years later Ehlers j^.40j gave transformations to derive 

such fields from vacuum ones. These are the only other papers on the subject written in non-relativistic units. Similar 
transformations were given by Bonnor ^ and with the help of the TWS method [Uliil. The latter was apphed to 
scalar fields and to the RN solution "plli^ . 

When there is no space symmetry present, the simplest harmonic function in cartesian coordinates is 



U (x,y,z) = l + y2^^^, (x-Xif + {y-yif + (z- Zif . (60) 

i 

It can be shown that the sources are point monopoles with masses mi and charges Si jl3l l2fl| | , connected by Eq. 
(54), which is true also in the Newtonian theory. It ensures the equilibrium between the electric and gravitational 
forces among the sources. Such multi-black-hole solutions satisfy certain uniqueness theorems [isl 11^ and possess 
axially-symmetric reduction with the points staying on the z-axis 47]. They have been generalized to expanding 
cosmological solutions [ZM . l49l | and the moving charged masses do not radiate electromagnetic or gravitational waves 

Solutions based on Eq. (60) are electrovacuum solutions and have singularities at the point sources. Increasing 
their number to infinity we expect to obtain charged dust in equilibrium. However, for this purpose one must add to 
the r.h.s. of Eq. (3) T^" — ij.c^u'^u'^ and introduce the current in Eq. (5). In this way charged dust with mass density 
fi and charge density a is described. In the general static case one can impose the Newtonian equilibrium condition 

±cr = \/G^, (61) 

which is the density analogue of Eq. (54), but here it holds for the sources of the field, not for a test-particle. Then 
one gets [sj 

/=(C + 0)^ (62) 

where C is some constant. Unfortunately, charged dust clouds are interior solutions and we cannot use asymptotic 
flatness to set C = 1. Neither can we turn off the electric field putting a = for then /i = 0. Neutral dust cannot be 
in static equilibrium - it collapses. Therefore, C is usually set to zero by a gauge transformation and then there are no 
root gravity terms. The situation is not clear, but the perfect square in Eq. (62) and the possibility of conformastatic 
MP solutions speak in favour of i? = 2 and consequently, D ~ k = in Eq. (38). Like in the multi-black- hole case 
one puts (j) — U^^, f — U~^, but U is not harmonic; it satisfies the non-linear equation [Tlll5ll | 

AU^±^aU^ (63) 

In the axially-symmetric case the equilibrium of charged dust was investigated exhaustively by Bonnor [s^ ] both in 
the Newtonian theory and in general relativity. Eqs. (31) acquire terms with 11 and a and serve as their expressions 
through (j) and /. The equilibrium condition reads 



1/2 



v^AtVv/7+(TV(/) = 



(64) 
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and leads directly to an arbitrary functional dependence j^l"^ = F {(j)). Inserting it into the field equations yields the 
analogue of Eq. (34) 

F {Fl - 1) + {FF^^ - Fl + 1) F^V(I)V(I) = 0. (65) 
When (f> is found, / is also known, while 

^^le'^ = FF^A(b + {FF^^ - ~ l) VcpVcj), (66) 

l^^tre^'^ = -FAcI) + 2F^W<j>W<j> (67) 

determine n and a. The analogues of Eqs. (32,38) become 

fc, = r - 1) ^^-2 (02 _ ^ k^^2r{F^^l)F-^^r(t>z. (68) 

One solution of Eq. (65) is F^ — 1. It leads to the MP solutions with dust, Eqs. (61-63). What is more interesting 
is that when F^ ^ 1 one can define 

^ = J {F^-lf'F-'d^ (69) 

and Eq. (65) becomes = 0, so 5* is analogous to but the function F (0) is arbitrary. After it is fixed we can 
find 4" and express everything in terms of it. Such solutions do not satisfy Eq. (61) and are singular [H^ . Some of 
them have, in addition, negative mass density, but in the remarkable example |52| of pure root gravity f — Cicj) there 
is a region of positive 

Charged dust clouds o f sp herical or spheroidal shape and obeying Eq. (63) have been studied extensively in an 
astrophysical context [54 . ISSl ISalSTllSq . They have a number of interesting properties when compared to usual stars: 
their mass and radius may be arbitrary, very large redshifts are attainable, their exteriors can be made arbitrarily 
near to the exterior of extreme charged black holes. In the spherical case the average density can be arbitrarily large, 
while for any given mass the surface area can be arbitrarily small. When their radius shrinks to zero, many of their 
characteristics remain finite and non-trivial. 

Some more general solutions of Eq. (63) have been given too [s^ |63|- They include the case of constant /i with 
U g iven in terms of a Jacobi elliptic function. The idea that may be concentrated on surfaces was also discussed 
|61|. Thin dust shells of spherical, cylindrical or plane shape were given as examples. Finally, a magnetostatic dust 
solution with functional dependence between / and A^p is also known '62]. 

The charge density required to satisfy Eq. ( 61) is quite small. It is sufficient that in a sphere of neutral hydrogen 
one atom in about 10^* had lost its electron [SJI- However, all charged dust models have one essential flaw: the 
equilibrium is very delicate and unstable and a slight change in cr would cause the cloud to expand or contract. On 
the other hand, the models discussed in the present paper do not depend on equilibrium conditions. They require 
just the creation of strong enough electric or magnetic fields. Of course, if the charged cloud is divided into parts 
with positive and negative particles (or ions), this can also induce EM-fields and root gravity. 

The appearance of the WMP relation (14) was studied also when the dust is pressurized, i.e., in the case of charged 
perfect fiuids [6^ l6^ and a gravitational model for the electron was put forth 16^. Like the case of charged dust, 
many different relations between / and are possible, but only partial results in the spherical case were obtained 
[66l |. Spherical charged perfect fluid interior solutions are reviewed in Ref. |6^ . 

The analogy between electrovac and stationary fields is also worth being mentioned. When there are no EM-fields, 
the equations following from the metric (56) become 68^ ,6£| 

2Au = -e-*"VrJV^2, AO. ^ 4ViiV17, (70) 
= 2uaUb + ie-4-r!,r!b, (71) 

= '^e'^'^y/lSabch!"', hab = iOa,b - (^b,a- (72) 
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The comparison between Eqs. (70-72) and Eqs. (57,58) yields the Bonnor transformation l4ll l69l | between 
electrovac and stationary solutions 

f = fl ^ = in, k = Aks, (73) 

given here for axially-symmetric fields. In the general case Eqs. (58,71) coincide after scaling the coordinates, while 
preserving ^ab !6p|- Eq. (72) simplifies for axial symmetry, 



-bJr- 



(74) 



For magnetovacs the transformation involves either \ — iO. oi simply x — which is a relation between the true 
potentials. One obtains imaginary EM-fields which should be made real by a choice of the integration constants. 

When the Bonnor transformation is executed upon Weyl electrovacs, one finds a differential system for the gravit- 
omagnetic potential 

LOr = —iripz, = iripr- (75) 

Its solution is Lo = r(^r, Cz = *V' with being another harmonic function. Hence, ip must be imaginary and the case 
D < should be used. Eq. (41) transforms into 

/-I = cosh + ^ sinh (76) 

with B — ib, —D — + A. This is exactly the Papapetrou (P) stationary solution [toI Flf. which has puzzled the 
researchers for ten years (1953-1963) with its absent mass termjtill the properly behaving Kerr solution was found. 

Stationary fields have been studied much more profoundly %|| because of their astrophysical importance. With 
the help of the Bonnor transformation many facts about them and their Ernst potential Eg = + ifl hold also for 
electro (magneto) vacs . For example, one of the degeneracies in R^^j^ leads to the functional dependence {lo) and 
the P-sohition [t^, ItsL ItI . [t^ , while another one leads to algebraically special fields [T^ . The class of metrics with 
shearing geodesic eigenrays and spin coefficient r = again gives the P-solution [t^ ■ There are several theorems Q , 
Sec. 18.7, about the subclass of conformastationary spacetimes III, 113 

ds^ = /, {dx° + uoadx^f - K'^ (x, y, z) {dx^ + dy^ + dz^) . (77) 

According to them, conformastats are of Petrov type D. All such metrics, however, are explicitly known and are 
axially-symmetric. On the other side, axisymmetric conformastats cannot have the so-called Ernst coordinates [80l | 
Eg and E* because Es = Eg {E*). This dependence is fs {ui) in disguise, so all such solutions belong to the P-class. 

Finally, it is known that when fs and O are not functionally dependent, Eq. (71) is more important than Eg. (70) , 
which follows from it |8lj . One can write a system of differential equations for ^ab and expressions for /, |82l l83j . 
both for a timelike or a spacelike Killing vector. In the axially-symmetric case an equation of fourth differential order 
for k results |8^. When applied to electrovacs, this shows that root gravity and usual gravity are complementary. In 
non-Weyl fields k (which is always proportional to k) plays the role of a master potential, like tjj (which is proportional 
to y/n) in Weyl fields. One is tempted to speculate that the purely electrically induced Weyl fields should bear no 
reference to k and the spacetime must be conformastatic {B = 2). 



V. ELECTROSTATIC EXAMPLES 



In this section four simple experimental set-ups are given where the effects of root gravity may be detected and 
measured. They involve plane, spherical, spheroidal and cylindrical symmetry. In these cases Weyl fields represent 
the most general electrovac solution and there is no ambiguity as to their appearance. 



A. The moving capacitor 

The study of plane-symmetric electrovac metrics goes back to 1926 [s^, but for a long time their Weyl nature 
remained unrecognized. A plane-symmetric example of a WMP field was give n first by Papapetrou 12] . Later Bonnor 
studied fields with (p ~ (/){z), which include also some non-Weyl solutions [13. Let us discuss the gravitational field of 
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a uniformly charged plane with charge density cr. The classic Maxwell potential (which becomes the master-potential) 
has both plane and mirror symmetry 

■iP = -2TTa\z\=q\z\. (78) 

It vanishes at the plane z = and goes to infinity when \z\ oo. In the case D = we should replace Eq. (78) into 
Eq. (39). Kar has proposed a coordinate system where <j) becomes harmonic, 4> — q\z'\. It gives a constant electric 
field and further deepens the analogy with classical electrostatics. Then the interval becomes 

ds^ = f {dx'^f - {dr^ + r^d^^) - [dz'f , (79) 

where / = (1 + q\z'\) . When D^Q one sees from Eq. (38) that k depends on r and in fact 

k = -f 9'^'- (80) 

The metric does not inherit the symmetry of its source. Now the Kar's gauge gives (p (z) — —B/2+q \z'\+a. Requiring 
flatness at z' = one obtains a — B/2 and Eqs. (40,41) yield 

ds^ = / - e^'i'-'^'y {r'dr^ + f-^dz'^) - r^r^d^\ (81) 

where f = I + 2qa \z'\ + q^z''^. The non-inheritance is obvious again. The D — Q case is obtained when a = ±1 
and is the only one with plane symmetry. This is our main argument that B = 2. Other values of B (including 
-6 = 0) introduce, in addition, the parameter a, which does not originate from the electric field, unlike the charge 
parameter q. Probably it is invoked by mass sources. As for the condition / (z = 0) = 1, its necessity will become 
clear in the following. In conclusion, the only plane-symmetric Weyl metric is given by Eq. (79). Another example of 
non-inheritance was given in Ref. |86|. where the metric is not rotationally invariant. The opposite situation is also 
possible - plane-symmetric metric, associated with semi-plane-symmetric EM-fields,87j. 

In order to obtain a regular global solution, one must satisfy the junction conditions |88| at the plane. The metric 
is continuous there. The extrinsic curvature reads 



Kaa = l\9z^\ ''\gaa),. (82) 



where a = x'^,(p,r and there is no summation. Its eventual jump at the plane determines the required energy- 
momentum tensor Sat of the massive surface layer 

^S%^rb-S%7: ^a^Kllt (83) 

a 

For the Weyl form of the axially-symmetric metric one has 

«;5°o = e""''(2w^-fc.)|t, K^^^^-e^-'^fc.lt, nS\^Q. (84) 

At first sight, the jump in the acceleration requires the introduction of mass on the charged plane. However, 
its true cause is the electric field present in the space around the plane and induced by the charge distribution on it. 
This is seen after one writes the Einstein and Maxwell equations, Eq. (31), in the present case 

u.. =e-2-02^ (j),,=2u,(j), (85) 

and makes use of Eqs. (14,37) 

u, = {l + <t>)^.- (86) 

Obviously, the jump in at z = is due to the jump in the master-potential tpz and consequently to the presence 
of charge. When we integrate the first equation in Eq. (85) across the thickness of the plane 5, the l.h.s. has a jump, 
no matter how small 5 is. The r.h.s. will be an integral of Tq which in the limit (5^0 coincides with S'g. Therefore, 
it also has a jump due to the charge. The functional dependence u (0) intermingles the equations in Eq. (85) and 
the charged surface layer of the plane causes the jump due to u in the extrinsic curvature too. For Weyl fields the 
usual procedure of introducing mass and pressures on the surface represents an attempt to model the influence of the 
electric field upon the metric by traditional mass sources [s^. l35l l89l|. When ■0 from Eq. (78) is replaced in gi, the 
plane will be attractive for one sign of 9, as if there was positive mass on it, but for the other sign of q it will be 
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repulsive as if it were made of exotic negative mass, which breaks the energy conditions. There is no paradox because 
the true creator of these effects is the energy-momentum tensor of the electric field, which always satisfies all three 
energy conditions. As a result, the term Uz in S% should be absent. In our case fcz = and even k — so there is 
no mass surface layer at all. 

Now let us put a second plane at z = d, charged in the opposite way. The electric field is confined between the two 
planes and V-' reads 

V' = ^^, E, = -fi:,. (87) 
Here -02 is the potential of the second plane (■01 = 0) . It is related to the charge density by 

1P2 = 2na-, (88) 

where e is the dielectric constant. Up to now we have considered the case e = 1. More generally, e enters Eq. (3) 
because the energy Too ^ ^Ef and consequently T^,y £Tfj_^. When absorbs the constants in Eq. (12) it will pick 
also ^/e. The same is true for ip. The acceleration formula (49) becomes in a dielectric medium 

9z = VgIJ^ « 2.58 X 10-^^i;2- (89) 
a a 

Taking finite disks instead of planes one obtains the usual capacitor. The field around its centre will be plane- 
symmetric while at the rim it will depend on r too. This may be diminished by careful electric shielding of the 
capacitor. Outside there will be a vacuum fiat metric which joins smoothly the interior due to the conditions / ~ const 
on the plates. One comes to the conclusion that the capacitor will be subjected to practically constant gravitational 
force Fg in the z-direction. 



Fg = VGi^V2 = \fG'etxS^2. (90) 



where M is the mass of the dielectric, (i is its mass density and S is the area of the plate. This force is very different 
from the electric force, trying to bring the plates together 

The latter is neutralized by the mechanical construction of the capacitor. If it is hanging freely, the effect of Fg may be 
tested experimentally. To increase the acceleration it is advantageous to make d small (typically 0.1cm < d < 1cm), 
to raise the potential difference ip2 between the plates up to 2 x 10'*CG5' and to take a material with high e. Some 
examples are gases (e w 1), quartz (4.5), glycerine (56.2), water, muscles, electric ceramics (81), rutile {Ti02) with 
e = 170. Ferroelectrics are even better if one can cope with their hysteresis and tendency for saturation in strong 
fields. Barium titanate (BaTiOs), potassium dihydrogen phosphate(i4ri?2f O4) and many others have e in the range 
of 10*. Thus yfefd may reach in principle 10^ and the maximum acceleration max = 5.2(7e is more than enough 
to counter Earth's gravity. In a more modest attempt one can take "02 = lOO/cF and \feld — 10^ to get about one 
percent of ge- Curiously, the same factors e,d,rp2 are even more important in Fe, which is always much bigger than 

As we shall see, root gravity effect is strongest in this first capacitor example and, more generally, when the electric 
(magnetic) field lines are parallel. The gravitational acceleration is not strictly constant, having an unobservable 
dependence on z. One cannot obtain root gravity terms by studying fields with constant acceleration |90l |. Earth's 
field also shares these "shortcomings" of artificial gravity: the acceleration's directions are not parallel but meet at 
the planet's centre and its magnitude decreases with height. In the traditional approach to the charged plane its 
metric depends directly on t and z, bypassing the axially-symmetric step, see Q, p234. The Weyl nature is then 
completely obscured. Some plane metrics are induced by null EM-fields, like the special pp-wave given by Eq. (15.18) 
from Q or the Robinson- Trautman solution, given by Eq. (28.43) from the same reference. They are non-static 
because null EM-fields are incompatible with static metrics, Q, Theorem 18.4, |93|. We also do not discuss the non- 
null homogenous and uniquely conformally flat Bertotti-Robinson solution, Q Sec. 12. 3. All other plane-symmetric 
solutions with non-null EM fields have been found by Letelier and Tabensky '95|. The metric is either static or 
spatially homogenous. We are interested in the static branch. It has been reviewed in Ref. [9^, however, the electric 
field was not discussed and the ties with the Weyl fields remained unelucidated. Let us clarify this issue now. 
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The metric reads in cylindrical coordinates 

ds^ = Kc^dT^ - N + R^d<^'^) - PdZ^, (92) 
where K, N, P depend on Z and there is one relation between them. It allows to express K and P as functions of N 



K 



1 + /32Va^) , p=0(;3i+/32V]v) \ (93) 

2 

/3i = |^, /32 = l-/3i, (94) 

where -q and /i are constants related to the charge and mass density on the plane. For the electric field we get from 
Eq. (5) 



We have chosen the boundary conditions (j) (0) = 0, (0) = 1. It is easily seen that 

if = l + Bo^ + f^', A^=fl + — 0) , (96) 



Thus if, the analogue of /, obeys Eq. (14) and root gravity terms are present. When is turned off by 77 ^ 
(/3i — > 0), Bq does not stay constant. It increases to infinity instead, so that i3o0 remains finite. This triggers the 
"mass out of charge" mechanism described in Sec. 3 and we end with the vacuum solution for a massive plane with 
K = N~^/'^ j^^l- Hence, the mass is not entirely of electromagnetic origin and the parameter /i is independent from 
77 in general. Another limit is (32 0. This was explored first by McVittie 95] . Then (3i — 1, rj — ±2/i, Bq — ±2, 

is: = i/7v. 

One can absorb N in Eq. (92) by the change Z ^ N. However, it is better to determine N from the relation 
between K, P and N which specifies the coordinate system. The conformastatic spacetimes studied by Weyl |Q| , 
Papapetrou T^] and Bonnor 19] have N = P. Kar sU\ used two gauges, KP = 1 and KP = N"^ in his pioneering 
work. The second gives constant Ez, as seen from Eq. (95). McVittie .95] worked in the gauge KP = N. Patnaick 
[9^ attacked the problem with K = P, which is the Taub's gauge in the vacuum case The same gauge was 

utilized by Letelier and Tabensky |9^ . Unfortunately, the equations cannot be integrated explicitly in this gauge 
except for the McVittie limit. 

In order to make comparison between the axially-symmetric approach and the traditional plane-symmetric approach 
one should use the gauge N = P. This differential equation for N is easily solved and yields 



iV = P^(l-^Z+^Z^) , 1-^Z) , (98) 



2N 



2fi 



l^^Z+^Z' 



^1/2 



- 1 



(99) 



These expressions coincide with the B ~ 2 case, Eqs. (39,78) only in the McVittie limit where (32—0 and consequently 
q = ji = 77/2, rj = — 47rcr. In this limit the fields in Kar's gauge [sM l93l ] coincide with Eq. (81) for a = 1 after the 
identification Z = z' . The original McVittie's metric is obtained by passing to e''^ = qz' + 1 



{dxy 



-2qz" 



(dx^ + dy^^ — e 



-iqz" 



(dz"f 



(100) 



In conclusion, the purely electric effect upon the metric of a charged plane is given by the McVittie solution, while 
the deviation of (32 from zero (/i from 77/2) signals the presence of mass in addition to the charge. In this way one can 
increase (very inefficiently) Bq and the root gravity acceleration. 



15 



The solution in the general plane-symmetric case is also axially-symmetric. There must exist a transformation of 
it to the Weyl interval, Eq. (20). The latter is obtained from the general static element (see Ch.8, Sec.l) in 
isothermal form 

ds^ ^ goo {dx°f + g^^dip^ + gii (dx^)^ + {dx^f . (101) 

The interval in Eq. (92) is a particular case of Eq. (101) and, as expected, V KNR satisfies the two-dimensional 
Laplace equation. Therefore, one introduces the two-harmonic coordinate 

r{R,Z) = ^/KNR = R[l-(^Z^ (102) 

and finds its conjugate z (i?, Z) from the condition r + iz = _F (i? + iZ), F being an analytic function. The Cauchy- 
Riemann conditions lead to 

z = Z /33 (i?2 - Z^) , (103) 

/ = r^R-^ (1 - + Ai?') , = (r" + ^PIR^Y^ , (104) 

where (3^ — /32M/4 while R^ is determined from the quadratic algebraic equation 

The addition of mass through the "mass out of charge" mechanism leads to fc 7^ 0, depending on both r and z, while 
/ acquires r-dependence. The plane Z = is transformed into the elliptic paraboloid of rotation z = /Jar^. The 
non-electromagnetic mass changes the shape of the charge-carrying surface when axially-symmetric coordinates are 
used. 

B. Repulsive gravity 

Let us study now the gravitational effect of electric fields with spherical symmetry. The master potential is 

= -, =r^ + z^. (106) 
P 

This is the exterior solution for a metal conductive sphere of radius pi, charged to a p otent ial V'l = q/pi- Inside the 
sphere = and the spacetinie is flat. Outside we obtain the charged j9^|93 Curzon |lO0j | solution. Eq. (38) gives 

^ = -^. (107) 

Like in the plane-symmetric case the metric does not inherit the spherical symmetry of ip unless D = (B = 2). This 
corresponds to the critically charged Curzon metric. Eq. (39) yields 

p-q V P P J 

The charge in CGS units q is connected to g by g = "^P?- Eq. (49) gives for the acceleration 

ffp = -yG7%l. (109) 

p 

It has a maximum at the sphere and utilizing our maximum potential we obtain l^p^maxl = 5.06/pi. When pi = 10cm, 
one gets 0.5cm/s^. Formula (82) for the extrinsic curvature holds in the present case after the replacement z —> p 
and a = x^ , 9, if, i.e., the axially-symmetric element is written in spherical coordinates 

ds^ = e^" {dx°Y - [e2fe (dp^ + p^de^) + p" sin^ 9d^^] . (110) 
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Expression (84) for 5 q and S"^ holds after the same change, while S'g = 0. Once again we argue that the Up term 
is in fact absent, while fc = and there is no mass surface layer, but only a charged one. The acceleration can be 
enhanced by introducing a second sphere with p2 > pi and filling the space between the two spheres with high e 
dielectric. A spherical condenser is obtained and the master potential is increased by y/ep2 {p2 ~ Pi) . As in the 
previous case, for one sign of q gravity becomes repulsive. The interaction betwee n massive bodies with repulsive 
(negative mass) and attractive (positive mass) gravitation has been discussed in Refs |lOll Il02j | . In our case repulsion 
is a natural property of the electric field and does not break the energy conditions. 

Let us deform the charged sphere into an oblate or prolate spheroid. Its gravitational field is described best in the 
corresponding spheroidal coordinates x, y (to be distinguished from the cartesian coordinates in the previous sections) 

r^T{x^±lf\l-y^Y'\ z^Txy, (111) 
where t is a parameter. In the prolate case one has 

X^^{l+ + l^)^-, y^^{l+-l_)^ (112) 



;± = ^(z±r)'+r2. (113) 

The master potential is taken to depend only on a; (a; xi is the surface of the charged spheroid). The harmonic 
solution is given by the Legendre function Qq (x) 

^^^In^^ (114) 

^ 2 x+l ^ ' 

and generates the general solution of the Einstein-Maxwell equations for such symmetry. From Eq. (48) one has 
q = Ana (x^ — l). This potential coincides in form with the Weyl rod for the vacuum 7-metric Eq. (38) gives 

Dq'^ , - 1 , , 

k=^\n- -, (115) 

8 x-^ — 

which depends on y too, signalling non-inheritance, except when B = 2. The same conclusion follows in the oblate 
case, but there ip = — g arctan l/x. Repulsive gravity also arises in these cases. 

It is well-known that the unique sp herically symmetric electrovac solution with mass M and charge Q is the 
Reissner-Nordstrom solution 

ds^ ^ + c^dT^ -(l-^ + 9l) \r^-R^ (de^ + sin2 ed*^) , (116) 

V R R I \ R R- 1 



GM 



= — , TO = 



GQ2 



Q2 = (117) 

Surprisingly, all three 
^ ^ ^ , criticall y (e xtremely) 

The proof uses the Weyl form of the RN solution HI ■ The 
essential step is to apply Kar's gauge, in which (j) becomes harmonic instead of ■(/;. When _D = we transform Eq. 
(108) by setting R = p — q, cos = z/p and obtain Eqs. (116,117) with to = — q, Q = q. The sphere p — pi transforms 
into the sphere R = Ri = pi — q. When Z) > wc take Eq. (114) and fix q by the condition Dq^ /A = 1. Then we 
utiHze Eqs. (42,43) to find and / 



Hence, at least the solution given by Eq. (108) should be transformable into the RN solution, 
solutions described above are formally equivalent to its three cases; undercharged {Q^ < ni^ 
charged (Q^ = w^) and overcharged (Q^ > to^). 



Making the identifications 



D \ ^ + x 



f = 



TO 



Q'f\ B 



x'-l 



2to 



2 • 



(118) 



(119) 
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we obtain (together with Eqs. (111-115)) the formulas for the undercharged case from Ref. The transformation 

R^Tx + m, cosQ = y (120) 

maps this solution into the RN solution. The charged spheroid x ~ xi goes into the charged sphere R~Ri~ rxi+m. 
It is seen from Eq. (119) that B oo when Q —^ and B(j) remains finite. The independent mass parameter m again 
arises from the "mass out of charge" mechanism. A similar chain of arguments connects the oblate spheroid case to 
the overcharged RN solution. 

In conclusion, the electrically induced spherically symmetric gravitational field is given by the critically charged 
Curzon solution, which is equivalent to an extreme RN solution with m = —Q. There is a transformation of RN 
solutions with m ^ ±Q into spheroidal metrics with particular D. The distance between the spheroid's foci r is related 
to the deviation of m from its electromagnetic value. Thus part of the exterior solution's mass is from electromagnetic 
origin. This has been long known for interior charged perfect fluid solutions ,^7] . The results are ana logous to the 
plane-symmetric ones. The RN solution also exerts a repulsive force for certain values of its parameters |l05lll06j . It 
appears, however, due to another mechanism. In the charged Curzon solution the sign of the charge is decisive and 
the region with repulsion occupies the whole exterior space. In the RN solution the mass is always positive, but enters 
the metric with a negative sign, so that a competition with the charge term is possible. The region with repulsion is 
finite, R < Q'^/m. 



C. The microgravity chamber 

It is natural to ask whether artificial gravity of electromagnetic origin can be created in a big enough chamber for 
people to work inside. This can be achieved by taking a capacitor and increasing the distance d between the plates to, 
let us say, 2m. The dielectric should be some gas, e ~ I and we cannot use its enhancing property. Since air breaks at 
Ecrit = 3 X W^V/cm = lOOCGS, we better use vacuum. Taking one creates microgravity, = 2.58 x 10~^cm,/ s^. 
One can compete with when ■02 = 3.8 x lO^'f/'niax- In addition, the exact solution is very complicated. Although 
the effect is tiny, it is worth to present a simple exact solution with the same properties. 

The general solution of the Laplace equation in spherical coordinates can be written, separating the variables, as 

V' = 5Z (^1"/^" + C2nP"^"+'^) Pn (COS0) . (121) 
n=0 

Let us take a charged sphere of radius pi . When n = 1 the interior solution has almost constant 

Tp~=qpcose. (122) 

The potential will be continuous when n — 1 also for the exterior solution 

+ qpl cose 3q 

tp^ = , a~ — cost*. (123) 

47r 

Here a is given by the jump of V'p at the sphere. Obviously ■0+ is a dipole solution. In classical electrostatics this 
charge distribution arises in a number of cases: 

1) Two oppositely charged balls with slightly displaced centres. 

2) Uniformly polarized ball. 

3) Metal uncharged ball in a uniform constant electric field. 

4) Electric dipole in a spherical shell cut inside a conductor. 

These set-ups are inapplicable for our purposes. One must create the cosine law density by using plastic or other 
non-conductive material, since the sphere is not an equipotcntial surface. Inside it is the plane-symmetric solution 
with B = 2 and k = Q. Outside k should vanish too, in order to be continuous. The gravitational accelerations read 

9z = vGq, g^ =0, 5j = gj = . (124) 

Outside the sphere there is a peculiar gravitational force ~ p^^ ■ The monopole term is absent and one may think 
that the source is a dipole of positive and negative mass particles. In fact, the source is the charge of the sphere and, 
as in the capacitor example, gravitation is concentrated mainly inside and rapidly decreases with the outside distance. 
Intuitively, this is a more economic way than to generate long-distance monopole terms. 
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The traditional junction conditions yield 

=2eXI-=2e2"Vp|t, S'g = = 0. (125) 

Again, the introduction of mass on the surface is not necessary - it will only duplicate the effect of the charge. 

Gravitational solutions corr esponding to electric or magnetic dipoles plus eventual mass distributions have been 
considered in the literature [sl llOTl ri08| . The global solution described here was discussed too ,109j . However, the 
spherical layer was made from dust, which requires for equilibrium the introduction of mass according to the law (54). 
Several definitions of acceleration were discussed in relativistic units, where the root gravity effect is obscured. The 
focus was put on the similarity between the interior solution and "gravity" in an accelerated coordinate system, and 
the radiation of test-particles. 



D. Charged cylinder 

Finally, we discuss the gravitational field of an infinitely long charged cylindrical shell of radius ri . In the exterior 

i^^gln-, V^. = (126) 
ro r 

where q = —Airari is determined by the charge per unit length a and tq is a scale not set by a. Inside the cylinder 
■0 = 0. The second gravitational potential reads 



4r 



e^^^r— . (127) 



The symmetry of ip is inherited because /, k depend only on r. The extrinsic curvature is obtained from Eq. (82) by 
making the change z ^ r. The same is true for 5"° and Eq. (84). Eq. (31) becomes 



^ -2u,2 ' 1 



— Ur — e "0„, 6rr H 0r — 2Ur(t>r- (128) 

r r 

Making use of Eqs. (14,37) we find an equation similar to Eq. (86), which shows that u^. has a jump induced by ipr- 
Tg also acquires a jump when Eq. (128) is integrated across the layer and it is exactly the Ur term in Sg. Hence, 
this term is redundant and the surface layer is given in terms of kr 

kS°o = kS'^^ = -e^-'^krlt, k5; = 0. (129) 

Put differently, a mass surface layer of this type (with equal mass and tangential pressure) generates a non-zero k in 
the metric, D < for positive mass. 

Studying purely electromagnetic effects on the metric, we set D = 0. We also demand that when the electric field 
is turned off {q = 0), we should have (j) = and / = 1. Therefore, in Eq. (39) i/jq = —1 and 

(j)= -1+ (^1-qln^^ , f = (^l-qln^^ . (130) 

The interval reads 

ds^ = (^l-qln—^ c^dt^ - (^1- q In— ^ {dz^ + dr^ + dip^) . (131) 

Comparing it to Eq. (22.16) from Ref. one sees that the latter refers to Eq. (22.4b) with A = 0, instead of 
Eq. (22.4a). Thus, t and z should be interchanged. Second, the a in Ai (corresponding to our q) should be in the 
denominator. Third, in Ref. 0] the condition '0g = was used and the gauge transformation 0' = 0-)- 1 was made. In 
this case / = 6'"^ a nd there is no root gravity term. However, when g ^ both the metric and 4> become ill-defined. 
Finally, Ref. IllOl is quoted for this result. In fact, Raychaudhuri applied the Rainich formalism to rederive some 
results of Bonnor [13 (who used the Weyl formahsm) and whose Eq. (2.30) he quoted. 

In reality the charged cylinder is also massive, so that one must allow for a tiny k to appear. Eqs. (127,129) show 
that a mass density rric induces negative D 

(132 
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It confirms our assertion that D is generated by "parasitic" mass sources. They can be neglected since kr n and 
the mass induced acceleration for realistic to^ is much smaller than the electromagnetically induced one. 

Let us discuss further the cases with D 7^ in order to make connections with the existing literature. We introduce 
the notation 



]-^\q, ^ = nln-, r± 
2 ro 



When > we use Eqs. (42-44) with -00 = to find 

^r, „ D 



B 
"2 



/ = 



\D\ 



= r: 



- ± - 



2 r_' ' 4sinh^<^ 
Using Eqs. (14, 37) one can express through / and 'tp to obtain in the present case 



B 
2" 



When I? < 0, Eqs. (40-41) with Vo = give 

B 1 

<b = \ — 

^ 2 2 



Dtan^ = 

2 



/^r_ {in) 
2 r+ {in) ' 



(133) 



(134) 



(135) 



(136) 



f^-—^^[T+{in)]-\ e2'==r-^"^ (137) 
4 cos^ (; 

These trigonometric expressions were given by Bonnor |19| and quoted by Raychaudhuri |llO| . We have explained 
already that the D > case does not simply follow from the substitution in n, as stated in Ref. [3 and Eq. (134) 
is an illustration of this. Eqs. (134,135) appear as Eq. (22.14) in Ref. where 



r~" r" B 

a, = -^, = Ao = <P+- (138) 



and b = —q, m = n. The credit is given again to Ref. jll(T| |. The same formula was re derived in a recent paper 
^llj. The particular case n = 1/2 was studied in the pioneering work on the subject |112| (see also Ref. |ll3j where 
the Rainich formalism was applied). In Ref. Q 'T' is considered real and the Bonnor's case D < is not mentioned. 
Yet, when a Bonnor transformation is made, it maps this case on the Lewis class of stationary cylindrical vacuum 
solutions, while the real n case {D > 0) is mapped onto the Weyl class (see Eq. (22.6) from Ref. 01). These two 
clas ses a re preserved also in the electrovac case. In addition, the metric (131) parallels the limiting Van Stockum class 
[7lllll4| . Like the D = case, the metrics (134,137) are singular because ^0 — 0- They become regular when i/jq is 
chosen according to Eqs. (40,44). The parameter ro in Eq. (131) can be set to tq = ri, making the metric continuous 
at the shell. Like in the previous examples, the cylinder's gravity is repulsive for one sign of a. 

It is easy to see that at r = the above solutions are singular, no matter how one chooses ro and D. This is the 
reason to use them only as exterior solutions. One can invoke the "mass out of charge" mechanism to end with the 
vacuum cylindrical Levi-Civita (LC) solution |^l23l| 



ds^ = r^'^c'dt' - r^"^ {dr' + dz') - r'~^-'dip'. (139) 

For this purpose we set VZ? = and send B and Z? > to infinity, keeping n finite. In this way one of the terms 
in T± vanishes and Eq. (139) results with n — 2v. On the other side, finite n means g — > 0, i.e. the electric field 
vanishes. B does not stay constant, consequently the term B(j) is preserved and from Eqs. (7,135) we get 

-^Ur^{-+<j>Ur^-^-^-. (140) 



The LC metric describes the gravitational field of an infinite line mass (ILM) and is singular at the axis. Eq. (140) 
shows that we have reproduced its acceleration out of a charged massless solution. 
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VI. MAGNETOSTATIC EXAMPLES 



We have pointed out that magnetic fields produce the same effects as electric ones. One must replace 4> by A, Ei 
by Hi, while tp becomes the ma gnet ic scalar potential, determined by surface currents. This is confirmed also by the 
definition of the magnetic field |115) 

= -\v^e,kiF''\ (141) 

which gives, using Eq. (22) 

= -A, = -fi,,. (142) 

Eqs. (21,37) provide the connection between the scalar potential and the only component (x = A^p) of the vector 
potential 

Xz = r'tpr, Xr = -r4)z- (143) 

It should be noted that in fiat spacetime magnetostatics the formula H = curia involves the physical component 
in curvilinear coordinates a^^-) = A^p/r. In magnetogravity one has to find the scalar potential anyway, in order to 
obtain /. 

In a general medium the energy is Too ~ I^H^, where /i is the magnetic constant and consequently T^i, — > ^T^^. 
Hence, picks also the multiplier ^/JI. Eq. (142) remains unchanged. Arguments, analogous to those for electric 
fields, lead to the conclusion that B = 2 for a pure magnetic effect upon gravity. Formula (49) for the acceleration 
becomes with high degree of precision 



g,^~^GfiH,. (144) 

It is similar to the electric case, Eqs. (53,89). The gravitational force is very different from the Lorentz force, acting 
upon charged particles 



Ft. = eE 



^(vxliy (145) 



In it the electric force does not involve e, while the magnetic does involve /x, but acts only on moving charges. There 
are terms, depending on the velocity, also in gi. We are interested on its effect upon macroscopic bodies for which 
f <C c so we have neglected them and study only acceleration at rest. 



A. Current loop 

The first magnetovac Weyl solution was given by Papapetrou . Later Bonnor presented a number of examples 
|21| . including the gravitational field of a current loop. He found a very small effect, based on mass-energy 
considerations. The root gravity term remained unnoticed. 

The magnetic field of a current loop is well-known and includes elliptic integrals. On the 2;-axis and at the centre 
the field simplifies 

77,(z,0) = 0.27r/ ^— -, Ho = H,{0,0)^- , (146) 

{rf + z^f'^ ri 

where ri is the loop's radius in cm, H is measured in Gauss and the current / is measured in Amps. Setting /i = 1 
we get from Eq. (144) 

30 = 15.(0,0)1 = 1.62x10-4-. (147) 

n 

Earth's acceleration is reached when He — 3.8 x lO^G = 380T or le/ri — 6.05 x 10^ A/cm. For a laboratory set-up 
let us take ri = 100cm. If a lightning with cross-section of radius ro = 10cm circles around the loop, one may take 
I = 10^ A, the current density being J = 3.18 x 10'^ A/cm^ and go = 0.162cto/s^. One can reach ge with a current of 
6.05 X lO^A. 
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It is advantageous to make the loop t hicker, tu rninK it into a finite solenoid. Let the inner radius be ri, the outer 
radius be r2 and the height be I. Then |ll6llll'<| 



Q, -(_ /q,2 _|_ 02 

Ho = F{a,P)riJ, F (a, /3) = 0.47r/31n ^ , ^ , (148) 

1 + vi + /^ 

where a = rijr\^ (3 — l/ri and J is the current density. As an example, let us take ri = lOOcm, r2 — 2ri, I — 2ri. 
Then « 1 and J = Hq/IOO. Now ge is reached when Je = 3.8 x lO'^A/cm'^. 

One can increase the acceleration by creating magnetic fields in a ferromagnetic medium. Iron has /imax = 5000 
(V^max = 70.7). There are alloys, like supermalloy (79% Ni, 16% Fe, 5% Mo), which have //max = 8 x 10^, y^/Imax — 
894.4. Their saturation field is comparatively low, Hg = 8 x lO'^G and the maximum is obtained roughly for one 
third of this value. The effective field in Eq. (144) will be Heff = -^Aimax^max ~ 238T, which is of the order of He- 
Making a disc from this material with a current flowing in a strip around the rim, we get the magnetic analogue of 
the moving capacitor example from Sec. 5.1. 

There are two conceptual issues which need clarification. The metric depends directly on the scalar potential, but 
ip is multivalued in magnetostatics. If there is a current-carrying surface, the jump takes place there and probably 
the metric can be made continuous in this region of non-vanishing J. In the case of a current loop the jump can be 
arranged to take place on any surface, based on the loop. Symmetry considerations require to take the disk z — 0, 
p < ri as such surface. In spherical coordinates |23| 



ip = 0.27rJ 



±1 _ iLft , lift _ „^ , y^3!L^ f^)""" p,„, , 

n 2rf 2.4.. .2n \ri J 



(149) 



where the sign coincides with signz and P„ {cos 9) are the Legendre polynomials. Now, let us remember that we have 
been unable to determine the sign of B and for definiteness worked with positive B. In fact, Eq. (39) should read 

f=il±^P)-\ (150) 

Using different signs for z positive or negative, and taking into account that P2n+i (cos ^) = 0, makes / continuous 
at z = 0. The derivatives of ip are single-valued, hence, k is also continuous. The acceleration will have a jump and 
change of sign at the disk within the loop, like it has on both sides of a charged plane (disk). We took the absolute 
value of z in Eq.(78), which plays the same role. This jump has been attributed to some mass surface layer |2lj . which 
in our view is not correct. 

The second issue concerns the fact that the magnetic field far away from the loop has a dipole character and the 
acceleration does not contain a monopole term. As we have argued, this is not a tragedy, because electromagnetic 
fields can create artificial gravity in a confined region of space. Intuitively, this is more favourable energetically. Under 
a Bonnor transformation such fields go into the Papapetrou stationary solution, which possesses the same feature. Of 
course, the material of the loop always has some mass, which induces a monopole term, negligible with respect to the 
strong effect from root gravity. 

B. Spherical solenoid 

The spherical solenoid is an example of a closed current-carrying surface and a magnetic analogue of the microgravity 
chamber discussed in Sec. 5. 3. The scalar potential is similar to Eqs. (122,123) 

^-=-Hopcoa9, = cos 61. (151) 

2p-^ 

Hq is the approximately constant magnetic field inside the solenoid, directed along the z-axis. The normal component 
tjjp is continuous at p ~ pi but ^jg has a jump according to Eq. (55) 

, I 1 - , I 47r , , 

Het = M- = —J- 152 

Pi c 



Eq. (151) gives 



Stt 

J = Jo sin 6', Hq = -^Joi (153) 
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where Jq is in A/ cm. This kind of current can be obtained by rotating a uniformly charged metal sphere or making 
a spherical solenoid with constant number of coils per unit length of the 2;-axis. The last condition ensures uniform 
magnetic field also inside a prolate or oblate spheroid. It can be fulfilled much easier than the cosine law for a in Eq. 
(123). Putting Eq. (153) into Eq. (144) yields 

= 2.17 X 10-Vo. (154) 

Earth's acceleration is reached when Jo = 4.51 x 10^ A/cm. This is better than a loop and worse than a fat solenoid 
because the sphere is thin, but the gravitational field is uniform and mimics the field of Earth. Outside the sphere 
the fields have dipole character and quickly vanish. The junction conditions again lead to Eq. (125) but now S'g = 0. 
No mass surface layer is necessary. 

This example was considered in Ref. j^. There the physical component was taken in Eq. (152), which is 
wrong, although the difference is negligible. The current determines the jump in the classical scalar potential, which 
is related to the classical magnetic field Hq. Perturbation theory was used in this reference, which starts with terms 
^ K and misses the root gravity term ~ -^k in the exact solution. It is no wonder that the results of Ref. |^ confirm 
the findings for a long solenoid jj, given here in Eq. (9). 



C. Long solenoid 

In order to explain why Eq. (9) is so different from Eq.(17), some preparatory work is necessary. Let us rewrite 
Eq. (31) in terms of /, A 

Aln/ = 2/-i(A2 + A2), AX = rHfrK + fzXz). (155) 

Next, let us introduce h = j j and express A through The second equation is satisfied trivially. The first one and 
the compatibility condition X^z — X^r give 

Alnh = -2h-^X^,+xl), Ax ^ h-' (hrXr + h^Xz) ■ (156) 

These equations differ from Eq. (155) only by a minus sign. Let us search for Weyl solutions in the system h, x by 
demanding h ~ h (x). The analogues of Eqs. (14,34) are 

f-7^-W Ax= ^ f A xl + xl). (157) 

Ci + Bix -X Ci+ Bix - X 

Assuming the fimctional dependence x (V') and Aip — 0, the analogues of Eqs (36,37) follow 

Xi = yV'ii = rip^, X^ = -rtpr- (159) 

The function / should remain positive when the magnetic field is turned off, hence, Ci > 0. It can be set to one, as 
we shall see. The analogue of D, Di = Bf + 4Ci is always positive, leaving only the option Di > for the integral 
in Eq. (158). This time, however, {Bi — 2%) < Di and one finds instead of Eqs. (42-43) 

Bi ^fWi , ^flh , ^ 4r2 2 



X = — ^ H ^ tanh ^-V', f = cosh^ ^li- (160) 

^ 2 2 2 £>! 2 ^ ' 

Since cosh a; > 1 always, / ^ oo when r cx) and Vo is useless. We put tAo = 0. Obviously, / does not fit for 
an exterior solution except, possibly, in the cylindrical case. When t/; = ■)/; (z) we have / (r, z) and already / breaks 
inheritance. For a regular interior solution at the origin, should compensate the factor r^. Using Eq. (32) with 
replaced by A, one obtains after some calculations 

fc=:ln-+fc, (161) 
r 

= - ^^l) , h = ^ri^r^z- (162) 
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Here k satisfies the same equation (38) as k in the Weyl case. Acceleration is obtained from Eqs. (7,157) 

9^ 5\ , x-Si/2 



r 1 + Bix — 



Xr (163) 

and obviously there is a root gravity term ~ B\\J Gxij"^- One can also express x in terms of / and 'i/'i from Eqs. 
(157,159) 

This formula is analogous to Eq. (135) in the electric case. 

When the magnetic field is turned off, the line element becomes 

ds^ = {dx"y - {dr^ + dz^ + dip^) . (165) 

According to one interpretation r and Lp become cartesian coordinates like z and this is the Rindler solution represent- 
ing flat spacetime in accelerating coordinate system. This means that the static source will be accelerating when we 
pass to a non-moving coordinate system and is not reasonable physically. The second interpretation is to consider Eq. 
(165) as the v = 1/2 case of the LC solution (139). There are well-known problems with the physical interpretation 
of this met ric w hen v > 1/2. Recently a source in the form of a cylindrical shell of anisotropic fluid was found for 
< z/ < oo |ll8j . Out of the several definitions of mass density per unit length, the most reasonable seems to be the 
Israel's one, which gives 

^^^"^ = A 2 — TT- (1^^) 

The mass reaches a maximum exactly for v = 1/2, this being a turning point. The coordinates z and Lp switch their 
meaning for v > 1/2. On the other side, the LC solution is flat for = 0, 1/2, oo. The general impression is that an 
infinite line mass field has been imposed upon a usual Weyl solution, endowing it with some cylindric features. Such 
metric with parasitic masses and singularities cannot compete with the genuine Weyl solutions except for one case. 
Let us find the magnetic analogue of the charged cylinder. One has 

V^ = i^oln-, ^r = —, (167) 
ro r 



f = r'rl, fc^n^lnr, n ^ "^f^ (168) 
^=TT^ + ^' H,=rA = Ho (169) 

tPrl + 2 

and Hr — 0. Eqs. (142,159) have been used, while r_|_ is taken from Eq. (133) with D Di. The physical magnetic 
field is 

12 / 

This solution is similar to the one given by Eqs. (133-135) and has the line element 

ds^ - r^T\ {dx'f - r^'^'rl {dr^ + dz^) - T^^dip\ (171) 

It is regular at r = and may be used as an interior solution only when n = 1. Then we put Tq = Di to obtain flat 
spacetime when Hq = and find 

ds2 ^ f [(dx°)^ - dr^ - dz^] - f-^r^difi^, (172) 

+ (173, 
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The electric analogue of this solution was discovered by Bonnor in 1953 1^, while studying solutions with 4> = 4> {z). 
Its magnetic version was found a year later |^. Later it was obtained implicitly in Ref. (119j . Melvin discussed 
its properties again in connection with Wheel er's t heory of geons. It is known as Melvin's magnetic universe, since it 
can be used globally. Gautreau and Hoffman |l2Clj gave the more general solution in Eqs. (160-162) and showed that 
the Bonnor solutions follow from it. It was also derived through coordinate modelling Here we have proved that 
this is the most general Weyl solution for the system (156). It is not a Weyl solution for the original system (155), 
but still has a root gravity term. The general cylindrically symmetric solution with Hz field is given in Ref. 1^ by Eq. 
(22.11). However, this is not in the axially-symmetric form because grr 7^ <7zz- That is why it differs from Eq. (171), 
although Eq. (22.13) coincides with Eqs. (172- 173) . A c ombination of LC and cylindrical electr omag netic solutions 
has been also studied in the Rainich formalism |l2lL Il22j and in the usual Einstein formalism 0, llllj . 

One can invoke the "mass out of charge" mechanism exactly as in the case of the charged cylinder. We set \/Di = 
and send Di to infinity, keeping n finite. The charge q is replaced by ifo- The LC solution (139) is obtained with 
n = 1 — 2v. This time Bix is preserved and one gets an acceleration formula from Eqs. (163,169) 

%^--^^ = -. (174) 
c'^ r Zr r 

The final result is the same as in Eq. (140), namely the acceleration in a LC gravitational field. 

After this digression, let us return to the exceptional case, Eq. (173). There are no root gravity terms in this 
formula and one may think th at at last this is the example of vanishing Bi. However, the regularity condition n = 1 
leads to Bi — 2\/ jG — H^/ Hq^ which, in general, is not zero. It is worth to trace, using the formulae from this 
section, how the acceleration loses its root gravity term, due to the regularity condition 

Bi\ 1 „ , 2rK.H^ 

^-T)^^^r+(1^^^+)^^32^T7^- ^'''^ 

We have mentioned that the Weyl cylindrical solution is always singular at the axis, hence, for an infinitely long 
solenoid the Bonnor-Melvin solution should be used. However, this does not seem to describe any physical situation. 
In any real situation, starting from the current loop, going through finite solenoids with open ends and ending with 
the closed-surface spheroidal ones, there exist a regular Weyl solution (depending on z or z and r) with root gravity 
term. The infinite solenoid appears to be the unphysical limit in a series of elongated prolate spheroidal solenoids, 
where 5z ^ changes into and root gravity into usual gravity. 

It should be noted that there is another solution for tp except Eq. (167), which makes / in Eq. (160) regular. 
Separating the variables in the Laplace equation in cylindrical coordinates, the radial dependence of ip is given by the 
Bessel functions Jq (r),A^o ('')■ The latter possesses a logarithmic singularity, which can compensate the term in /. 

The transfer of the system (155) into the system (156) has, due to the Bonnor transformation, its mirror among 
statio nary vacu um solutions. It is given by Theorem 19.3 from Ref. based on Ref. 123] and rediscovered several 
times |l24l ll2,'Tj . The functional dependence in Eqs. (157,158) translates into the so-called S (A) solutions (0, Sec. 
(20.4) 

S = r\-^^^ ^u;^ + B2u; + C2, i^s ^ J (176) 

where Alps = 0. They are also known as the Ehlers class 126"! and were investigated recently |127| . The elucidation 
of this connection comes to show that the study of S (A) solutions began already in 1953. 



D. Line current 



In Maxwell magnetostatics when a current / flows in a straight line conductor along the axis z with radius ri, it 
creates a circular magnetic field in the outside 

H^ = —, (177) 

cr 

where CCS units are used, the current included. It follows either from a scalar potential ip {(p) or from a vector 
potential with only a component 

21 21 

tp = (fi, ttz — Inr. (178) 
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In four-dimensional notation Oz = Ar.- In general relativity, instead of defined by Eq. (141) one should use in this 
case the physical component i?(<^) |l3j| 

ii{ip) = P(z){r) = (179) 

In flat spacetime it yields the usual definition H = curia. Such type of field breaks the circularity condition, necessary 
for the usual projection formalism of stationary (and static) axisy mmetric fields (Q, Sec. 19. 2). Nevertheless, in the 
cylindrically symmetric case the solution is known explicitly [l28l| (see also Eq. (22.11)) 

ds^ = r2"'r2 (^{dx°f - dr^) - r^VH^^ - T'^dz^ (180) 



I ^ f F(^^)(^,.)y/gvv9ttd(pdt (183) 



e = r = air" + a2r-", 03 = 40102^^ > 0, (181) 

where ^ = A^. The metric of the infinitely long solenoid (171) also can be written in this form when the changes 
Az A:p, z ^ ip are made. Eq. (180) is wrongly quoted in Ref. |^ as referring to the long solenoid - the Melvin 
metr ic, Eq. (22.13), does not follow from it. One should consult the first edition o f Rcf. 3] on this question. After 
|l28j| . the gravitational field of a line current was studied also in Refs. |llllll2lLll22j . Combining Eqs. (179,180), one 
obtains 

H(^) = (182) 

In order to generahze Eq. (177) one should use Eq. (10.87) from Ref. Ts^ (see also |l22j ') 

47r 
c 

Here < t < 1 and the circular integral holds for any r > ri. This gives 

^r-r^'^H^^)-^- (184) 

The flat spacetime limit follows when n = 0, F = 1. The field is non-Weyl, but has a root gravity term. Indeed, we 
can use the analogy with the long solenoid's equations (157,168) to obtain 

r-^ = l + B,C-e, 900 = ^^^^^^^2 - (185) 
Then Eq. (7) gives the analogue of Eq. (163) 

|.2!,(,_f)r»,. ass, 

Replacing Eq. (184) into this equation gives for the dominating root gravity term 

gr = -O.IBiVG-, (187) 
r 

where / is in Amps. The metric (180) is singular at the axis, but anyway we use it as an exterior. Therefore, the 
regularity mechanism leading to the Melvin magnetic universe is not needed here. Unfortunately, there are no hints 
for the scale of Bi. Assuming its typical value 2, finally yields 

gr = -5.16x10-"" — . (188) 
n 

For a "static" lightning / = lO^A, ri = 10cm, gi = 0.516cm/s^. The magnitude of ge is obtained for, let us say, 
ri = 1cm, le = 1.9 X lO^A. The gravitational force is radial, like in the long solenoid example, does not follow the 
magnetic lines, but is perpendicular to them. 
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VII. CONCLUSIONS 

When the derivations, proofs and other technicahties in this paper are omitted, the summary of the results remains, 
divided into 19 points. Some of them are not new and the corresponding references are cited in these cases. 

1) The gravitational acceleration at rest in Weyl-Majumdar-Papapetrou fields has a root gravity term, proportional 
to = -s/G, which is 10^^ times bigger than the usual perturbative coefficient (?k. It is linear in the EM-fields, 
while the perturbative term is quadratic. Sizeable gravitational force exists (see Eq. (53)) although the metric is very 
close to the flat one (Maxwellian limit). Its explicit form determines up to a sign the important constant B. For its 
typical value B = 2 the Earth's acceleration is obtained in electric fields of order lO^V/ cm and in magnetic fields 
of about 380T. One can change the direction of gi by changing the direction of Ei or Hi and confine gi to a finite 
volume by confining the EM-fields. 

2) In WMP fields the gravitational potential depends directly on the four-potential of the EM-fields. This is 
not amazing - such potentials are known to be important in certain quantum effects. One should not make gauge 
transformations because they spoil the asymptotic behaviour of the metric or the initial condition - when the EM-field 
is turned off, flat spacetime should result. Coordinate transformations are formally equivalent to modification of the 
symmetry of the charge (current) carrying surface and to a change in the physical set-up. 

3) The energy-momentum tensor (in particular, its energy component Too) induces a change in the Ricci tensor ^ k 
according to the Einstein equations (3). This leads to changes in the metric and its acceleration, which can be ^ ^/k, 
and may contain no monopole term. Creating artificial gravity that is localized in space and has no long-distance 
mass terms seems energetically more favourable. 

4) In axially-symmetric systems Weyl fields provide regular exterior and interior solutions to any distribution of 
charges or currents on a closed surface. They are determined by a master-potential '0 satisfying the Laplace equation. 
When the metric depends on just one coordinate (three commuting Killing vectors are present) the Weyl solution 
becomes the most general one. In truly axisymmetric cases it presumably determines the pure electromagnetic effect 
on gravity, while the solutions of the Ernst equation include hidden mass sources. The constant B (taken positive 
for definiteness) divides the Weyl solutions into three classes, according to B = 2, B > 2, B < 2. 19]. Among them 
S = 2 is privileged, being the simplest (conformastatic spacetimes). When B 2, parasitic masses appear, either on 
the surface through the junction formalism, or on the axis. In some cases the metric does not inherit the symmetry 
of the EM-source. 

5) The gravitational force at rest induced by electric and magnetic fields is the same 

[Hill 113 

unlike the Lorentz 

force, acting upon charged particles. The surface sources determine the master potential and not 0, . 

6) There is a "mass out of charge" mechanism, which allows to obtain solutions with mass and charge from the 
Weyl solutions. It clearly indicates the part of mass which is of electromagnetic origin. Eq. (14) still holds and root 
gravity term remains, but B is affected by the mass. Such solutions can incorporate the mass of the charged surface, 
which is always present in practice. 

7) In the general static case a harmonic master potential for WMP fields appears only when B = 2 H^ . 
Point or line sources have singularities, hence, one should use closed surface sources (shells). Another alternative 
is to use volumes of charged dust or perfect fiuid, where the functional dependence / (0) appears naturally as an 
equilibrium condition, but is arbitrary [53 . For charged dust a harmonic potential can always be introduced |52| |. but 
the equilibrium is unstable |63l l64l l65l l6a | . It is more realistic to charge conductive surfaces by applying potentials 
or pass currents through coils wound around them, relying on the motion of free electrons in metals. One can also 
create powerful fields by separating the positive and negative charged ions. 

8) The wealth of solutions to the stationary vacuum Ernst equation is preserved completely for the static elec- 
tro(magneto)vacs due to the Bonnor transformation (41j,|6^]. Usual Weyl solutions correspond to the Papapetrou 
stationary solutions, while Weyl solutions to the transformed system of equations (156) correspond to the S {A) so- 
lutions. They consist of only one class D > Q and although non-Weyl, possess root gravity terms. The case S = 2 is 
special also for stationary solutions. 

9) In order to construct global solutions around the charged shells, the fulfilment of the junction conditions is 
required. They have some subtleties in the Weyl case. Eq. (14) interrelates the Einstein and the Maxwell equations 
(31) and a jump in T^^ results due to the jump in tp, caused by the charges or the currents. The term in Sab, 
containing Ui is absent and a mass surface layer appears only when ki is non-trivial. The acceleration is caused by 
the EM-field and not by often unrealistic fluid sources, trying to duplicate this effect. 

10) The pure electric plane-symmetric effect on the metric is described by the McVittie solution, which is a Weyl 
solution of class B = 2. When B ^ 2 the Weyl metric is not plane-symmetric (non- inheritance). Solutions with 
mass and charge also contain a root gravity term. There is a coordinate transformation between them and a Weyl 
axisymmetric field, which changes the plane into a paraboloid. 

11) A usual freely hanging capacitor is an excellent tool for testing root gravity. Constant acceleration is induced in 
its dielectric, Eq. (89). Taking maximum potential difference, minimum distance between the plates and a material 
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with very big dielectric constant (ferroelectric) one can obtain = 5.2(7e. This force is trying to move the capacitor 
in the z direction. It is smaller than the force of electric attraction between the plates. 

12) The pure electric spherically-symmetric effect on the metric is described by the charged Curzon solution, Eq. 
(108), which is related by a simple transformation to the critically (extremely) charged RN solution. This is a 
Weyl solution with B — 2. When B ^ 2 non-inheritance follows. The same is true for prolate and oblate spheroidal 
solutions. They are singular at the axis, so one can use them as exteriors to a charged spheroidal shell with Minkowski 
interior. Their gravity becomes repulsive in the whole outside region for one sign of the charge. Typically, the repulsive 
g — O.bcm/s". There is a coordinate transformation between RN fields with mass and charge and Weyl fields outside 
spheroidal charged shells. 

13) An exact Weyl solution exists with practically constant gravitational force inside a charged sphere, which quickly 
decreases outside, according to a dipole law. Unfortunately, with potentials mankind can create at present, the force 
inside this microgravity chamber is very weak, although still detectable. 

14) The pure electric cylindrically-symmetric effect on the metric is described by Eq. (131), which is regular outside 
the shell. In non-regular form it was first given by Bonnor |0| and is the electric counterpart of the limiting Van 
Stockum stationary solution. Once again, B = 2. When B 2, a mass surface layer appears on the shell, inducing 
the second gravitational potential k. The "mass out of charge" mechanism leads to the vacuum LC solution. 

15) The magnetic field of a current loop induces gravitational force, Eq. (147). The Earth's acceleration is reached 
when the current / = 6 x 10^^. A thick loop (fat solenoid) requires for this goal current density of 3.8 x W^A/crn^. 
Presumably, a ferromagnetic disc with a current strip around it is the magnetic analogue of the moving capacitor. 
The ambiguity in the sign of B seems to be necessary to make the metric unique and continuous across the plane of 
the loop. 

16) A spheroidal solenoid represents the magnetic analogue of the electric microgravity chamber. The uniform 
gravitational force inside equals that of Earth when J = 4.5 x 10^ A/cm. Outside the fields quickly vanish, being 
dipole in character. 

17) The Weyl magnetic solutions of the transformed system of equations (156) are magnetovac mirrors of the 
stationary S (A) solutions. They have singularities due to parasitic infinite line masses, consist of only one class and 
possess root gravity terms. In the cylindrically-symmetric case the metric is similar to the charged cylinder, but it is 
regular at the axis when a special condition is satisfied (Bonnor-Melvin solution, also used in the interior of a infinitely 
long solenoid f3|). The root gravity term is sacrificed in favour of regularity, resulting in an acceleration ~ kHq, which 
is negligible for realistic magnetic fields on Earth. The infinite solenoid seems to be unphysical limit, because in any 
finite solenoid there exists a regular Weyl solution with root gravity term. 

18) The cylindrical gravitational field of a line cur rent breaks the circularity condition of the axial symmetry 
formalism, but nevertheless can be found explicitly 0, Il28l | . It is similar to the field of a charged cylinder or to the 
field of the magnetic solution from the previous point. It is non-Weyl, but still has a root gravity term. If i? = 2 a 
current of magnitude 1.9 x 10^ A creates a radial gravitational force, equal to the Earth's one at r = 1cm. 

19) Some statements in the literature have been corrected. The Weyl case D > does not follow simply by analytic 
continuation from D < 0. A number of misunderstandings in the cylindrically-symmetric Einstein-Maxwell fields in 
Ref. has been clarified. 

VIII. DISCUSSION 

It is hardly believable that such an important property of WMP fields as root gravity has been overlooked in the 
past. If we consider the work of Reissner as implicitly pioneering, WMP solutions date back to 1916, the year of the 
Schwarzschild solution. In the past 88 years the subject has been explored pretty thoroughly as can be seen from the 
number of references. Only in the first half of the previous century root gravity could have been discovered 8 times: 
by Reissner, 1916 103], Weyl, 1917 [11, Nordstrom, 1918 104], Kar, 1926 ;85|, McVittie, 1929 [H, Mukherji, 1938 
\n^ . Papapetrou, 1947 12] and Majumdar, 1947 jll|. There were plenty of opportunities later too. In the same 
time certain spherically-symmetric perfect fluid solutions have been rediscovered up to 7 times. 

One reason is the wide-spread use of relativistic units - already Weyl worked in them. The only papers on the 
subject (except the present one) written in non-relativistic units are by Ehlers [s^li^- Unfortunately, they are in 
German and are extremely rarely cited. When G = c = 1, -^/k/Stt ~ k/Stt = 1 and there is no difference between 
root and usual gravity. When SnG = c = 1 then k = 1 and n/dtir ~ 0.2 while k/Stt = 0.04. The difference is just 
an order of magnitude. 

Another reason is that the easier solutions with plane, spherical, spheroidal or cylindrical symmetry were studied 
directly and were almost never treated as subcases of the axisymmetric solutions. This cut their ties to Weyl solutions 
from the very beginning. The RN solution was shown to belong to the Weyl class only in 1972 [4^ . 
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A third reason is the shaky position of the hnear term in Eq. (14), which can be ehminated in principle by a gauge 
transformation. In many papers disappearance has been its fate. We have argued that physicaUy this is not possible 
and the typical value of i? is 2. 

A fourth reason is the underestimation of exact solutions in favour of approximation schemes and numerical tech- 
niques. Root gravity escapes undetected by these methods. In fact, it sets a new scale ^/k or \fG and perturbations 
should be done around the exact Weyl solutions. 

More generally, WMP solutions have never been in the main trend of development of general relativity. Black holes, 
perfect fluid star models and cosmological solutions write its official history. The fact that the unique charged black 
hole is a WMP solution pushed their invest igati on towards multiple-charged-blac k hol es plL IT^ l20l liH 143 . WK liS 
H^H^i particle trajectories in these fields |l29l | or even WMP-based wormholes |l30l |. A similar astr ophysic al s tring 
rings in charged dust clouds (Bonnor stars) and their non-singular interpolation to black holes 0, EE 113, IH, Il3ll 
Il32j | . Axisymmetric solutions have been flooded with generation techniques for the Ernst equation. Papers on WMP 
solutions often appeared in local, old or unpopular journals, unavailable to the scientific community. Even the Weyl 
seminal paper still remains not translated into English. 

A sixth reason concerns the motivation of the first researchers. Weyl himself was much more interested in his 
conformal theory of gravity, trying to unify gravitation and electromagnetism. He never returned back to the Weyl 
fields. Ironically, McVittie presented his solution as an example to test one of the unified theories of Einstein. The 
father of relativity used to complain about the difhculties of exact solutions and found only four of them (Einstein 
static universe, cylindrical waves, spherical cloud of counter-rotating particles and the Einstein-Strauss vacuole). He 
concentrated his energy on the geometrical unification of the two fundamental long-range interactions. According to 
him, the l.h.s. of Eq. (3) was made of marble, while the r.h.s. was made of wood. His goal was to replace wood by 
marble. 

In this paper we share another view, closer to the Rainich " already unified theory" . We are interested in the 
electromagnetic effects upon gravity. It is not necessary to use his formalism, which i nclu des products of Ricci tensors 
and is similar in complexity to the Gauss-Bonnet terms in the string effective action [133]. If is wooden, let us do 
some alchemy, let us obtain marble from wood, gold from lead. The point is that electromagnetism has been mastered 
for centuries on Earth and we know how to create strong EM-fields. As for gravity, we only know how to watch this 
force in action around the Universe. 

Our approach is not to explain but to construct and parallels the efforts to build wormholes, time machines and 
warp drives. We use, however, fields which are common and satisfy all energy conditions and not exotic matter or the 
Kasimir effect. Second, we try to induce gravitational force without altering very much the metric. This is possible 
thanks to the 20 orders of magnitude collected in c^/2 in Eq. (17). We have demonstrated that except the Newtonian 
limit, general relativity possesses also a Maxwellian limit. Root gravity is one way to do this, but there are others 
too. The next step should be to change considerably the metric of spacetime and invoke typical Einsteinian effects. It 
seems more difficult and expensive, but perhaps among the 450 line elements given in Ref. 3] or in the 1822 references 
quoted there, there waits the necessary solution. 

Our tool is well-known electromagnetism. It is hard to imagine that the strong or the weak interactions, which deal 
with quark confinement and particle decays, can be used to manipulate gravity. At least root terms are not possible. 
QCD and QFD interact through non-Abelian gauge fields and in the r.h.s. of Eq. (1) there will be cubic and quartic 
terms in A^. This makes impossible to swallow the Einstein constant as this was done in Eq. (12). 

It see ms th at the first who studied experimentally the interaction between gravity and electromagnetism was 
Faraday |l34j . He dropped or lifted heavy bodies inside solenoids or measured their charge after such movements. He 
expected that gravity would induce currents or charges, but the results were negative. His last paper, which described 
these futile experiments, was written in 1860 and rejected by Phil. Transactions. Two of the decisive experiments 
confirming general relativity included the effects of static masses on null EM-fields (light bending around the sun and 
the redshift of photons in a gravitational field). It is time to do the opposite. It is the electromagnetic effect upon 
gravity that is much more interesting. Unfortunately, the main device for strong magnetic fields - the usual solenoid, 
is often considered infinite for simplicity, which results in a solution without root gravity terms '3, T^l . We 
have given here several laboratory set-ups where root gravity can be detected and general relativity tested once more, 
this time in its highly non-linear part. The most promising seems to be the well-known capacitor and its magnetic 
analogue. The effect has not been discovered yet, because capacitors are used in low-voltage circuits and they are 
first firmly fixed and charged afterwards. Small root gravity effects should be present also in today's solenoids with 
strong magnetic fields. The acceleration is masked by that of the Earth, being a fraction of it. 
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